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The full quantum-statistical theory of the Vertical-Cavity Surface-Emitting Laser 
(VCSEL) in the form of the Langevin equations is constructed for arbitrary relations 
between the frequency parameters. The same theoretical treatment as in Ref. [1,2] 
are used. For detailed analysis the theory is applied for lasers with equally living laser 
levels and on this basis the analytical expressions for the spectral densities of the 
Stokes parameter fluctuations are obtained in the explicit dependence on the physical 
phenomena, including the spin-flip and the optical anisotropy. It is demonstrated the 
arbitrary distribution of electrons between the sub-levels under pumping does not 
restrict a possibility to achieve the noise reduction below the quantum limit. Under 
comparison with phenomenological treatment Ref. |3| it is shown this approach turns 
out to be not quite satisfied. 
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I. INTRODUCTION 

In the last years there has been an increased interest to the polarization properties of the 
VCSELs. This interest is motivated in the first line by the potential applications of this type 
of lasers in the high-rate optical communications Ref. Q|. But there is also more fundamental 
reason for understanding of the polarization behavior in VCSELs, namely, a possibility of 
generating the intensity-squeezed light using the sub-Poissonian pump of the active medium 
Ref. 0, |J To date, squeezing in VCSELs has been demonstrated experimentally for both 
the single-mode operation and in the multi-transverse-mode regime Ref. Q|. 

Now in the different scientific groups two models of laser are discussed in the main. First 
of them is with the shortly living upper level and the other - with the equaled lifetimes for 
both the levels. The known spin-flip theory |8( was elaborated for the first time just for the 
latter system. Nevertheless a lot of it is suitable for the first one too. 

As for statistical aspects the laser with the equally living levels was studied only phe- 
nomenologically Ref. |3f in distinguish from the other system which has been studied in 

n q 

details within the limits of the quantum electrodynamics Ref. |2|. Our main goal here to 
develop not less qualitative theory for the system with the equally living levels. This will 
allow us not only to write the different correct analytical expressions, but also to estimate 
the possibilities of the phenomenological treatment Ref. 

The paper is organized as follows. In Sec. |H]the basic equations of the theory will be given 
in the adiabatical approximation. In Sec. IIII| IIVI the fluctuations of the Stokes parameters 
will be introduced into consideration, be written the respective linearized equations and 
their solutions. In Sec. 0the comparison of our results with the phenomenological ones will 
be made. At last, in Sec. EU the possibilities for the polarization squeezing observation will 
be discussed. 
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II. THE LANGEVIN EQUATIONS FOR THE TWO-LEVEL LASERS WITH 

TWICE DEGENERATED LEVELS 



The well known spin-flip VCSEL theory Ref. [q] is based on the representation about the 
semiconductor medium as a two-level system with a twice degeneration of the levels (see 
fig. 1). One pair of sub- levels is connected by wave with the a + -polarization and the other 

- with the <t_ one. The upper "atomic" states \a±) relax with the rate 7 a , and the lower 

- with the rate 7&. For the description of overturns of the electron spins (the spin-flip) we 
introduce into consideration some incoherent processes \a + ,b + ) — > |a_,6_) and back with 
the respective rates ry( a > b \ 

We will not mention all the details, how the main equations are produced. They can 



be found in the previous work Ref. l], l2L l2(- The main items are as follows. At first, the 
equations for the Hermitian operators of the field amplitudes with the <r±-polarizations, the 
populations of the four sub-levels and the polarizations of the two actual transitions are 
constructed on the basis of the full quantum theory for field and matter. The obtained 
equations (the Langevin-Heizenberg ones), first, are operator and, second, contain the oper- 
ator sources (the inhomogeneous terms) in distinguish from the respective dynamical theory. 



Usua 
Ref. 



lv for writing the correct correlation functions one uses so-called Einstein relationships 



id- 



Next, to make the mathematical situation much simpler, one passes to c-number repre- 
sentation: all the normally ordered operators can be converted onto respective c-number 
9I lic|. As a result the equations, called the Langevin's equations, read: 



functions Ref. 

a± = — na± — (k s + uop) a T + gP±, (2.1) 

P± = - (7± + iA) P ± + g(N a± - N b± )a± + F P ± (t) , (2.2) 

N a ± = fi a ~ laN a ± - 7 < a) (iV a± - N aT ) - g(a* ± P ± + a ± P*) + F a± (t) , (2.3) 

Nb± = fi b - %N b ± - 7 W (jVfct - N bT ) + g(a* ± P ± + a ± P^) + F b ± (t) . (2.4) 

Here a± are the a± complex amplitudes. The atomic c-number variables are represented by 
the polarizations P± and the populations N a ± and N b ±. One can see the theory turns out 
to be very complicated mathematically, because generally speaking there is a system of the 
12 differential equations relative to the 12 variables. 

The frequency coefficients have the physical senses: k is the spectral width of the cavity 
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laser mode, g is the atom-field coupling constant, fi ajb is the mean rate of the incoherent 
pump to the upper, lower laser level, 7^ is the mentioned above constant of decay of the 
upper, lower level, the spin-flip rate, 7j_ is the rate of the transverse atomic relaxation, 

A = uo — v - detuning of the laser frequency v from the frequency of the laser transition u, 
the coefficients K a and uj p present the linear dichroizm and linear birefringence, connected 
with the optical anisotropy of the semiconductor crystal. 

As for the stochastic sources in the Langevin equations F b ±, F a ±, Fp ± , their properties 
are given by the following non-zero correlation functions: 

[ la + T W) N a± + /i a + ^ c a) N aT - g (a*±P± + a±Pt)- 



F a ±(t)F a± {t>) 

-PaHa/%] S(t 
F a ±(t)F aT (t>) 



(2.5) 
(2.6) 



F b± (t)F b± (t>) 
-p b fM b /2} 5{t-t' 



(76 + 7^ } ) N b± + in + ^N bT - g (a* ± P ± + a±P|)- 



F b± (t)F bT (t>) = - ji h \N b± + N bT )+ Pbf i b /2 6{t - t' 



(2.7) 
(2.8) 



F a± (t)F b± (t>) = g (a* ± P ± + a±P|) 5{t - t'), 



(2.9) 



F* P± (t)F P± (t') 
Fp±(t)F P ±(t') 



27± - 7a - 7c 



f)AT a± + /i a + 7 ^ ) iV (lT l 



■ 2g (a±P±) <y(t-f), 
(7 & + 7f) 1^)^-0 



(2.10) 
(2.11) 
(2.12) 



F P ±(t)F b ±(t>) 

Here the parameters p aj 6 determine a statistical aspect of the pump: the Poissonian (quite 
random) pump takes a place with p ab = 0, the sub-Poissonian (strictly regular) - with 
p a>b = 1. Under the calculation of these formulas it was proposed that distribution of the 
electrons between sub- levels turns out to be random even with p a b = 1. 

The Langevin equations (j2.1j) - ()2.4|) with the correlation functions (|2.5jl - (|2.12jl are suitable, 
generally speaking, for description of any two-level laser, because there is no any restriction 
on the relationship between the different frequency parameters. Further we will make the 
area of our interest narrower adopting our theory for the VCSELs with the equally living 
laser levels. 
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As mentioned in Sec. H] the two VCSEL models was developed. In one of them, which 
is studied in the San Miguel's group (including the pioneer work Ref. the lifetimes of 
both the laser levels were equaled. In the other model (the Giacobino's group Ref. 
to the contrary the lifetime of the lower level was much less than the lifetime of the upper 
one. As stated above we want to apply our theory for the first model and further we will 
put: 7 a = 7^ = 7 and 7^ a - ) = 7^ = 7 C . Besides it is suggested the incoherent pump takes a 
place only to the upper level: fi a = \i (fi b = 0). 

To compare our calculations with the ones, presented in Ref. Q , we introduce new func- 
tions D and d instead of the populations N a ± tb ±: 

D= l - \{N a+ + N a _) - (N b+ + N b _)] (2.13) 
d = \ [(N a+ - iV a _) - (N b+ - N b „)] (2.14) 

and put that the transverse relaxation constant is the highest between the others and we 
have a right to apply the adiabatical approximation. In the approximation the equations 
read: 

d± = —na± — (« a + iu p ) a T + c (1 — la) (D ± d) a± + £± (2-15) 
D = fi - jD - 2cD (| a+ | 2 + | a_ | 2 ) - led (| a + \ 2 - | a_ | 2 ) + ^ D (2.16) 
d = - ls d - 2cD (| a + | 2 - I a_ | 2 ) - 2cd (| a+ | 2 + | a_ | 2 ) + f d (2.17) 

Here 

2 

c= J- 2V a= V -^ 7s = 7 + 2 7c . (2.18) 
7 ± (l + a^) 7± 

The stochastic sources in the equations ()2.15|) - (j2.17|) are some linear combinations of the 
initial ones. The respective expressions can be found in Appendix [XJ 



III. FLUCTUATIONS OF THE STOKES PARAMETERS AND 
LINEARIZATION OF EQUATIONS 



The usual approach to non-linear equations is to try linearizing them relative to some 
small parameter (s). In the cases of statistical theories of laser systems these small parameters 
are introduced as additions in the exact solutions to the stationary solutions. 
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We will study here only one stationary regime of generation with the linearly polarized 
emission. Certainly, everywhere further we must imply the set of the physical parameters 
which is able to ensure a stability of this solution Ref. Ill I. 

nn 

Following to Ref. pi the stationary semi-classical solutions can be written in the form: 

a±st. = Q e~ iAxt , A x = uj p + aK Xl k x = k + n a (3.1) 

Here the value Q is expressed via the pump parameter r = /// n t h. {l^th. - the threshold pump 
rate) as 2Q 2 = r — 1. 

The equality between the circular components means we have the linearly polarized along 
the x-axes non-zero solution. Really, because the Cartesian field components are expressed 
via the circular ones in the form: 

a x = -^={a + + a_), a y = -^={a + - a_) (3.2) 

a x = V2Q, a y = 0. 

The respective stationary solutions for the active medium are: 



P±*. = -(l-t<x)D st .Q, D st . = JL— , d st . = (3.3) 
g 7 + AcQ 2 



According to (|2.16|) in the absence of the stochastic source the condition of the stationary 
lasing (D = 0) in the explicit form is: D st = k x /c. This provides us with a direct connection 
of the mean radiation power Q with the mean pump rate of the medium \i: 

— (3.4) 



7 + AcQ 2 c 

We have discussed the stationary solutions of the problem and now the exact solutions with 
taking all the fluctuations into account read: 

a± = (Q + 6a±(t)) e~ tA *\ D = D st + 8D(t), d = d st + 6d(t), (3.5) 

For the VCSELs we focus on the polarization effects. That is why it is convenient to 
introduce into consideration the Stokes parameters instead of the complex field amplitudes. 
They are expressed via the Cartesian field components as: 

So = K| 2 + Kl 2 ' S l = \a x \ 2 - \a y \ 2 , S 2 = a* x a y + a x a*, S 3 = i(a* x a y - a x a*) (3.6) 
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Respectively in our case of the linearly polarized regime their fluctuations are: 

5S 1 = V2Q(Sa* x + Sa x ), 8S 2 = V2Q(6a* y + 5a y ), 6S 3 = iV2Q(6a* y - Sa y ) (3.7) 

( 8 So — 5 Si) and the linearized system of equations reads: 

SS! = 7 (r - l)5D + (3.8) 

SD = -tyrSD - 2kJS 1 + £ D (3.9) 

and 

5S 2 = 2n a 5S 2 - 2u p 5S 3 - « 7 (r - 1)5 d + £s 2 (3.10) 

8S 3 = 2n a 5S 3 + 2u p 8S 2 - 7 (r - 1)<M + ^ 3 (3.11) 

<M = -T s 5d + 2kJS 3 + r s = 7s + 7 (r - 1) = 2 7c + 7 r (3. 12) 

One can see there are two independent systems of linear differential equations. The rela- 
tion between the stochastic sources ^ S4=123 , £d, £d and the initial ones can be found in 
Appendix |XJ 

IV. THE SPECTRAL DENSITIES OF FLUCTUATIONS OF THE STOKES 

PARAMETERS 

Introducing the Fourie image of the function F(t) as 

+oo +oo 

F n = -= f F{t) e tnt dt, F(t) = — f F n e~ lQt dSl (4.1) 
V27T J y2n J 

— oo — oo 

we can rewrite the equations (j3.8|) - (j3.12|) as algebraic ones: 

-tSl S3 lta = 7 (r - 1) 5D n + £ Sl . n (4.2) 

+ 7r ) 5Dq = _ 2 k x 5S liU + £ An (4.3) 

and 

- (zfi + 2n a ) 5S 2in = -2uj p 5S 3:Q - « 7 (r - 1) 5d + £ Sa>n (4.4) 

- (ifi + 2K a ) 5S 3 ,n = 2cu p <5S 2 , n - 7 (r -l)6d + £s 3l n (4.5) 

+ r s ) 5d n = 2k x 5S 3 , n + ^ n (4.6) 
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Now relatively simply to solve it, and the interesting solutions expressed via the stochastic 
sources in the explicit form read: 



5Sx 



l{r - 1) £,d,q + {ir - £ Sl ,n n 2 + itt^r - 2 7 ^(r - 1) 



-i 



(4.7) 



SS 2 ,n = 1/2 7(r - 1) (u p + an a + iVta/2) £ diQ - (u p (T s - zfi) + ajK x (r - 1)) £s 3 ,n+ 



+ (- (T s - til) (« a + iQ/2) + 7^,(r - 1)) £s 2 ,n 
(u 2 p + K 2 a + ittn a - fi 2 / 4 ) ( T s - zfl) - 7(r ~ ( K a - «cj p + ifl/2) 



(4.8) 



5S 



3,n 



1/2 -7(r - 1) (auj p - K a - zO/2) £ d>n + ^ p (T s - zft) £ &i 



- (r s - zQ) (« a + «fi/2) £<j 3l n x 

(cj 2 + + ifl Ka - fi 2 /4) (r s - zfi) - 7 (r - 1)k x (« a - au p + ifl/2) 



(4.9) 



So-called spectral densities of the Stokes parameter fluctuations (5 Si SSp.)^ take an important 
part under physical discussion. They are defined under writing the spectral correlation 
functions: 



SSifl SS k ,n' = {5Si 5S k )n S(fl + Q') 



(4.10) 



To get the wished spectral densities we must make the last preliminary step and write 
the respective spectral densities for the stochastic sources. According to Appendix |X] and 
(|3.1|) - (|3.4j) they are given as: 



(&)n = (&)n = £&)n = (& C5 3 )n = ^(r 2 - 1) 7 /c (4.11) 
(&)n = «xr(r-p/2) 7 /c (4.12) 
(d)n = ^r( 7s / 7 + r-l) 7/c (4.13) 

We define the non-zero spectral densities of the stochastic sources in the same way as for 
the Stokes parameters ()4.10|) : 



£,D,d,n£,D,d,n' = (£b,d)n + 



Si.n' 



(4.14) 
(4.15) 
(4.16) 
(4.17) 
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Now it is not difficult to get: 



(6Sfh = n xKx (2Q 2 (r + 1) + 7 2 r(p - pr + 4))/Ai, n 
(SS 2 ) n = 2n x K x ({l 4 /4: (r + 1) + a 2 {l 2 + b 2 )/\ n 
(5S 2 ) n = 2n x K x (n 4 /A (r + 1) + a 3 ft 2 + 6 3 )/A n 
W^n = 2n x K xl {r - l)(a 23 Q 2 /A + b 23 )/X n 



(4.18) 
(4.19) 
(4.20) 
(4.21) 



where the coefficients are defined as: 



a 2 = {u 2 + K 2 a + r 2 /4)(r + 1) + 7 (r - l)(a 2 riy4 + (au p - K x )(r + 1)) 
b 2 = 7 2 (cu 2 + n 2 a )(r + 1) + 7s7( r - 1) + « K a) 2 + 2n(au p - K a )(r + 1) 



+ 



+7 2 (r - l) 2 n 2 (r + l)(a 2 + 1) - (cj p + an a ) 



a 3 = (^ + n 2 a + 7 2 /4)(r + 1) + (aa; p - T s /2)(r - l)(r + 1) 7 + 7 (r - l)rT s /4 
6 3 = T 2 (r + 1)(k 2 + ^ 2 ) + 7 (r - l)r s [r(a 2 ^ 2 - K 2 a ) + 2K a (auj p - K a ) 



a 23 = a 



r s + 2(r + !)(«- K a 



b 2 3 = Is r(uj p + a>K a )(auj p - K a ) + (r + l)(u p (au p - K a ) + k(uj p + a« )J 
+7(r - 1) (u p + an a )(au] p - K a ) + (r + l)(a 2 + l)/ca> f 



+ 



Xi.n = fa 2 - 2 7 «:,(r - I))' + fi 2 7 V 2 



fi 2 



H 2 /4 + k 2 + cu 2 - Ka7s + 7 (r - 1)(k - /c a )/2 + fi 2 (/t a - T s /4)+ 



+7(r - !)«.,. (aWp - « a ) + T s (k 2 + w 2 



V. COMPARISON WITH RESULTS OF THE PHENOMENOLOGICAL 

CALCULATIONS 



As mentioned above the phenomenological calculations were made in Ref. j2| and the 
respective curves were drown as a result of some numerical analysis. Here we are going to 
compare some of these curves with our respective ones and then to estimate the phenomeno- 
logical approach on this basis. 
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In the fig. 2 and 4 are presented the spectral dependences, which we have copied simply 
from the cited Ref. |3j. The curves have drown with the following set of parameters: 7 = 
1GHz, k = 300GHz, u p = 1GHz, a = -3, 7, = 100GHz, r = 1,04 in absence of the 
linear dichroizm Kn = 0. 

Accord,,, Re, 13 the curves on F ig . 2 ar e the b „ dependences of the full spectra! 
power of the laser emission (the curve with single maximum) and - of the circularly polarized 
components (the curve with two maxima). These match to our (SSf)n and (5Sj)n + (5S%)q, 
which are constructed with help of the formulas ()4.18jl - (]4.19p and presented in fig. 3. 

As is seen our and obtained phenomenologically curves turn out to be qualitatively alike 
especially as the positions of maxima on the frequency axes are exactly the same. At the 
same time one can see the levels of noises in our approach turn out to be incomparably 
lower than in the phenomeno logical one (please, take into account, the vertical axes on all 
the pictures are chosen in the logarithm scale). 

The frequency dependences for the normalized cross-correlation 

(5Sf) n - (6Sp n 
° + - (5Sf) u + (5S!)u 1 ) 

are presented again respectively for phenomenological (fig. 4) and our approaches (fig. 5). 
Fixing here some qualitative (and quantitative - relative to the position of maxima) similar- 
ity, nevertheless one can see some serious differences take a place. For example, the right side 
of the graphics after maximum exhibits the absence of correlation in the phenomenological 
calculation and the appreciable anti-correlation about —1/2 according to our formulas. 

It is easy to understand why the positions of the maxima are the same along the x-axes 
in different approaches. It connects only with so-called relaxation oscillations and that is 
why this is perfectly independent of the choice of stochastic sources. 

As for quantitative differences we see the two reasons, why they take a place. First, 
under phenomenological introducing the sources to get the respective correlation functions 
one must use some additional physical considerations. Usually, especially under conditions 
of non-linear field-matter interaction, it is not a simple problem and requires a special 
attention. Regrettably this side of the question has quite fallen out of discussion in Ref. 0| , 
and certainly it provides us with a possibility for doubts relative to the obtained correlation 
functions. At any rate we have to fix our formulas for the correlation functions of the 
stochastic sources are quite different. 
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Second, one of the reasons of our differences is very clear. In Ref. Q the Langevin's 
equations were produced under introducing the sources into not the complete system of 
the dynamical equations but into the simpler system, where already two variables D and 
d appear instead of four ones N a ± and Nj,±. We remember use of the simpler system of 
equations turned out to be possible in the dynamical theory of the VCSELs with the equaled 
relaxation constants of both the levels. But in the statistical theory it leads to the losses of 
the important sources, which depend on the populations of the sub-levels N a ± and N b ±. We 
think it is main reason why the correlation functions in Ref. jjj are proportional only to the 
population differences of the kind of N a ± — iV&±, what can not be correct. 



VI. POLARIZATION SQUEEZING IN THE VCSEL EMISSION 

In Ref. J3| we have discussed in details how to observe the Stokes parameter fluctuations 
in experiments with two photodetectors. Choosing the available geometry of the experiment 
we are able to select the wished signal. The most important cases read: 

(5t 2 _) n /(t + ) = 1 + 2 K /n x (SSl) n , y = 0; 
(5ii) a /(i+) = 1 + 2K/n x (5Sl) n , f = tt/4, 9 = 0; 



{6i 2 _) a /(i+) = l + 2K/n x {6£fl) l 



(p = tt/4, 9 = tt/2; 



(6.1) 
(6.2) 
(6.3) 



and 



ip = n/4, 9 = tt/4 : 

(5z 2 _W(z+> = 1 + K/n x [(5S 2 2 ) n + (5S 2 3 ) Q + 2(5S 2 5S 3 ) n 
(p = 71-/8, 9 = 0: 

{5i 2 _)n/{i + ) = 1 + K/n x [(6Sl)a + (6S 2 2 ) n + 2(SSxSS 2 ) u 
<p = tt/8, 9 = tt/2 : 



(Si 2 -h/{i + ) 



K/n x {5Sl) n + (5S 2 3 ) n + 2(SS 1 SS- 



3)n 



(6.4) 



(6.5) 



(6.6) 



Here p is the angle between the direction of the linear polarization of the VCSEL emission 
and the polarization beam-splitter axis, 9 is additional phase shift, introducing by the phase 
plates between the orthogonal field components. 

Our theory elaborated in the previous sections gives the possibilities to study any signals. 
But here we will discuss only the first of them (|6.1J1 . connected with the polarization squeez- 
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ing and that is why having the principal character for quantum optics. As for the others to 
our mind it is interesting to consider them together with the respective experimental date 
for comparison. 

In semiconductor lasers (including the VCSELs) it is easy achieved the regularity in the 
pump, that, as we remember, leads to the essential intensity noise reduction below quantum 
limit for two level lasers without any degeneration. The degeneration of the laser levels 
introduces to system some additional random process, namely electrons under pumping are 
distributed between the sub-levels quite accidentally. This turned out to be inessential for 
the lasers with the shortly living lower level, because there in the linearly polarized emission 
only the full population of both the upper laser sub-levels N a+ + iV a _ plays the role for 
squeezing. And just it does not fluctuate under the regular pump. 

At the same time in our case the situation appears more complicated for understanding 
because now the value (N a+ + N a _ — N b+ — N b _) produces squeezing. 

To make some conclusion it is enough to watch only over the point Q = in the spectrum 
(jfi.lj) . It means we want to watch over a depth of the noise reduction below the quantum 
limit. Putting p = 1 (the regular pump) and n a = (without the dichroizm) we obtain that 

(« 2 _W<H> = 1 + '~ j^j ^ \ (6.7) 
Comparing it with the respective formula for the laser without the level degeneration j^j]: 

(<^W<H> = 1 + ~ ^ — \ (6-8) 



one can see there is an additional factor in 1)6.7)1 equaled to r/(r — 1). It is clear this plays 
an important role for a small amount of the pump parameter r (especially just near the 
threshold with r-l<l). At the same time, squeezing is able to appear only with the high 
enough pump parameter, namely with r > 5. Then the additional factor is already about one 
and hence turns out to be quite inessential. So our general conclusion is the degeneration of 
the laser levels does not leads to some additional difficulties in the production of squeezing. 

Also it is interesting to compare two kind of the VCSELs with the shortly living lower 
level (the case in Ref. j2|) and with the equally living levels. The respective formula in 
Ref. I2J is: 

(Si 2 - h=o/(i + ) = l + 3 -^-^-^0 (6.9) 
r — 1 r — 1 
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For example, if we choose r = 6, then according to the last formula we have the noise level 
is 0,28 of the quantum limit. At the same time for our case here it is only 0,84. One can 
see the reduction of noises below the quantum limit is more effective for the VCSEL with 
the shortly living lower laser level. 

VII. CONCLUSION 

To conclude we would like to say once more which concrete targets have been achieved 
in this work. First of all, the full quantum-statistical theory of the VCSEL with the equally 
living laser levels has been built. Thereby now we have a possibility to consider the problem 
on the same level of understanding as for the VCSEL with the shortly living lower level. 

On this base the spectral densities of the Stokes parameter fluctuations have been written 
under taking into account the optical anisotropy of the semi-conductor crystal (the linear 
dichroizm and the linear birefringence) and the spin-flip. 

We have been discussed the experimental situation in which it is possible to observe any 
spectral densities of the Stokes parameter fluctuations and also their correlations. Under 
discussion our main attention was devoted to the problem of polarization squeezing and 
the role of the degeneration of the laser levels. We concluded that the role of the random 
distribution of the electrons between the sub-levels under pumping is inessential for the shot 
noise reduction. 
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APPENDIX A: STOCHASTIC SOURCES IN THE LINEARIZED LANGEVIN 

EQUATIONS 



The basic equations of the theory ()2.1j) - (j2.4j) have the stochastic sources as inhomogeneous 
terms. The properties of them are specified by the non-zero correlation functions ()2.5j) - (J2.12j) . 
The case of our interest is [i a = fj,, /j,^ = and 7 a = j b = 7 and 7^ = = j c . We can 
rewrite the equations in the form which they is given in Ref. |3( in: 



a± = —kcl± — (n a + lujp) a T + gP±, 

P± = - (l± + iA) P ± + g(D ± d)a± + F P± (t) , 

D — 11 — jD — g(a* + P + + a+P; + a*_P_ + a_P*_) + F D (t) , 

d = -j s d - g (a* + P + + a+P; - alP_ - a_P*_) + F d (t) , ls = 7 + 2 7c 

Here instead of the populations N a ± and N b ± the new variables are introduced: 

1 



D = - 

2 

<-5' 



(iV a+ + iV a _) - (iV b+ + iV 6 _) 
(N a+ - N a _) - (N b+ - N b __) 



Respectively instead of the initial sources F a ± and F b ± the new ones are created: 

Fd = - (F a + + F a - — Fb+ — Pfe-) 
Fd= - {F a + — F a - — P&+ + Pb-) 
The correspondent non-zero correlation functions read: 

F D {t)F D {t') = i [ 7 (F a+ + F a _+¥ 6+ + iV b _) + 2//(l-p) 

-4# (<p7 + ^p; + ^!7 5 r + 5F )] 6(t-t'). 



F d (t)F d (f) 



(7 + 4 7c ) (iV a+ + iV a „ + iV b+ + iV b _) + 2/i- 
-4#(^^ + ^Pf + ^P7 + 0^)] 5{t-t'). 
F D {t)F d (t') = F d (t)F D (t>) = ~ [7 (iV a+ - iV a _ + iV fe+ - F 6 _ 
-4# (<P7 + ^Pf - Zl 5 ! - aZPTj] 5{t-t'). 



F D (t)F P± (t') = -- (7 + 7c) P± - 



(Al) 
(A2) 
(A3) 
(A4) 



(A5) 
(A6) 



(A7) 
(A8) 



F d (t)F P± (t') = T - (7 + 7c) P± S(t - t'). 



(A9) 

(A10) 

(All) 
(A12) 
(A13) 
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F P± (t)F P± (t>) = [(2 7i . - 7) N a± - 7c (N a± - N aT ) + fi] 5(t - *'), (A 14) 

F P± (t)F P± (t>) = 2g (^P±) S{t - f), (A15) 

Under the adiabatical approximation P± = 0, and we can get: 

1 



[g (D ±d)a± + F P ±] 



(A16) 



71 + 1A 

Taking this into account we can write our basic equations in the adiabatical approximation 
(|2.15|) - ([2.17|) . There the sources are: 



£ = - (1 - 1a) F P± 
9 



(A17) 

(A18) 
(A19) 

The linearization of the equations relative to the fluctuations near the stationary semi- 
classical solutions leads to the equations ()3.8j) - (j3.12j) which the fluctuations of the Stokes 
parameters are introduced in and correspondent ly the new sources read: 



£ 5l = V2Q{i x e^ 1 + £e 



-iA x t\ 



S 2 



■iA x t 



* a —iA x t\ 



£ Ss = tV2Q(C y e-^ - i y e 

£ D = F D -V2Q(£ X + C) 
Z d = F d -iV2Q(( y -Q 



iA x t\ 



(A20) 
(A21) 
(A22) 
(A23) 
(A24) 



where 



IV2 



iV2p 



;i - za) (F P+ - F P _ 



(A25) 
(A26) 
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FIG. 1: Four-level configuration of the VCSEL medium 




v [GHz] 

FIG. 2: The full spectral power (fragment of the figure 1 P|)with the following set of parameters: 
k = 300GHz, Up = 1GHz, ; a = -3, 7 = 1GHz, j s = 100GHz and n a = 0, r = 1.04 
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FIG. 3: The full spectral power (numerical calculation in our approach) with the following set of 
parameters: k = 300GHz, uj p = 1GHz, ; a = —3, 7 = 1GHz, 7^ = WOGHz and K a = 0, r = 1.04 
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FIG. 4: The frequency dependence of the correlations between the circularly polarized field com- 
ponents (fragment of the figure 1 with the following set of parameters: k = 300GHz, oj p = 
1GHz, ; a = -3, 7 = 1GHz, 7s = 100GHz and K a = 0, r = 1.04 
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FIG. 5: The frequency dependence of the correlations between the circularly polarized field 
components (numerical calculation in our approach) with the following set of parameters: k = 
300GHz, uj p = 1GHz, ; a = -3, 7 = 1GHz, 7s = 100GHz and n a = 0, r = 1.04 



